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ABSTRACT 


As global efforts advance toward a more environmentally sustainable future, there is a significant surge in 
the availability of electric vehicles for purchase. A growing number of electric vehicles does, nevertheless, 
give rise to novel challenges, including the requirement for environmentally sustainable and efficient 
charging alternatives. At this juncture, the microgrids initiate operation. Decentralized energy systems, 
such as microgrids, have the capability to function autonomously or in tandem with traditional centralized 
energy infrastructures. It possesses the ability to operate autonomously or in contrast to the conventional 
system. There are several benefits associated with microgrids, such as increased productivity, reliability, 
and resiliency. This article aims to evaluate the feasibility of integrating microgrids into electric vehicle 
charging stations. Case studies of successful microgrid deployments will be analyzed, the challenges and 
limitations of using microgrids in charging stations will be assessed, and the benefits of employing 
microgrids in this context will be discussed. Upon perusing this article's conclusions, you will possess an 
enhanced comprehension of the potential for microgrids to furnish environmentally sustainable and 
efficient charging solutions for electric vehicles. In Tangier, the scarcity of electric vehicle charging 
stations can be attributed to the substantial expenditure required for infrastructure development, such as 
electrical grids. On one side, the current infrastructure lacks sufficient development, thus impeding the 
widespread adoption of charging stations. On the other hand, Morocco has made significant progress in the 
field of renewable energy, particularly through investments in solar and wind power. Consequently, this 
study holds great potential for yielding valuable insights. 

Keywords: Microgrids, Renewable Energy, Energy Management Plan, Electric Vehicles, Charging 

Stations 

1. INTRODUCTION A significant number of the sizing studies for 
hybrid PV-WT systems have been created and 


A dependable and efficient charging infrastructure 
is becoming increasingly crucial as the prevalence 
of electric vehicles increases. A promising solution 
for powering electric vehicle charging stations is 
microgrids, which are decentralized power systems 
capable of operating autonomously or in 
conjunction with a larger grid. In light of the 
growing demand for electric vehicles in Morocco 
and the potential cost savings, resilience 
enhancements, and energy efficiency improvements 
associated with microgrid implementation across 
residential, commercial, and industrial sectors, in 
addition to charging stations, are all viable 
applications. In addition to their potential to 
revolutionize the way electric vehicle charging 
stations are powered, this essay will examine the 
application of microgrids in a the charging station. 


based on heuristic optimization approaches such as 
the Genetic Algorithm in [1], [2], [8], and [16], the 
Bat Algorithm in [11], and the hill climbing 
algorithm. In [4, 6], [9], [10], and [15], several 
surveys about the various ways of sizing were 
carried out. 

In this paper, the optimal scaling problem of a 
hybrid photovoltaic system used to power a 
charging station of electrical vehicle installation is 
investigated. MILP (Mixed Integer Linear 
Programming) is utilized in the formulation and 
resolution of such a problem. The problem of sizing 
has been solved using a variety of techniques, 
including the cutting plane method, Branch and 
Bound, Branch and Cut, Branch and Price, and, of 
course, the heuristics previously mentioned. For 
additional information, see [18] and [19]. 
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This study examines the charging station for 
electrical vehicles installation cost estimation based 
on MILP. Our methodology differs from that of 
previously published studies [17], [7], and [13] in 
that it incorporates fewer binary variables. Note 
that this reduction in the magnitude of binary 
variables is crucial for the proposed MILP 
problem's solvability. A simulation test comparing 
the Genetic Algorithm (GA) and the Branch and 
Cut Algorithm (BCA) is used to demonstrate the 
effectiveness of the proposed method. Our 
simulation demonstrates that BCA is superior to 
GA. 

The structure of the subsequent sections is as 
follows. A PV-wind turbine hybrid system is 
described in detail in Section 2. Our constrained 
sizing puzzle is described in Section 3. Section four 
of the paper presents the optimized hybrid model. 
Case studies are discussed in Section 5. Section 6 
concludes with concluding remarks. 


2. SYSTEM DESCRIPTION 


The system under investigation is comprised of 
Wind Turbines (WT), Photovoltaic cells (PV) and a 
battery energy storage system (BESS) are both 
components of a photovoltaic system. Figure 1 
illustrates this system's architecture. 

The solar energy captured by the PV cells is 
converted into electricity, which is then used to 
meet demand and/or charge the batteries. AC/DC 
converters are installed between the PV 
panels/BESS and the load because the energy 
generated by the PV panels is DC energy. This 
connection is facilitated by the Point of Common 
Coupling (PCC). 

2.1 PV System illustration 

The following formula is used to calculate the 
output power of a multi-crystalline PV system 
based on previous research published in [20] and 
[21]. 

P(t) = nAGr (t) (1) 


In which: 


e 1 is the total efficiency derived from n= 
nPV nSYS 
e Ais the PV panel's surface area. 


e GT is the total amount of radiation per 
square meter per hour. 
nSYS is the system efficiency, which depends on 
the efficiency of the AC/DC converter and cable 
losses, and P V is the PV module efficiency. 


The efficacy of the PV modules is dependent on 
temperature variations, as represented by the 
equation below: 


nev = nL! — B(Te — T:)] (2) 


Where nr is the rated efficacy of the module. 


e Bis the efficiency temperature coefficient. 


e Tr is the temperature at which efficacy is 
rated. 
Regarding the cell temperature Tc, it is computed 
based on the instantaneous ambient temperature Ta 


NOCT-—T, 
To(t) = ——#4 G(t) (3) 
GT.wacT 
With: 


e NOCT represents the Normal Operating 
Cell Temperature. 


e Gryocr is the solar radiation in NOCT 


2.2. Wind Turbin Model 

In previous works, the expression of the 
energy output of a Wind Turbine is regarded to be 
proportional to the cube of wind velocity (see, for 
example, [23] and [24] where this power is 
expressed). can be calculated using the formula 
below: 

Pw) = Ye Cp pAV %0) (4) 
where 
* p is the air density in kilograms per cubic meter. 
* A is the swept area of the rotor in m2. 
* Cp is the Wind Turbine's power coefficient, which 
describes its efficacy 
¢ V (t) represents the wind speed in m/s. 


2.3. BESS Model 

To control the transfer of energy from and 
to the BESS modules, it is essential to have 
information about the level of energy available in 
the batteries at each time t. This information is 
represented by the State Of Charge (SOC) 
parameter. 
The expression used to determine the state of 
charge can be found in a variety of sources, 
including [3] and [5]. These works demonstrate that 
the State of Charge variation can be calculated 
using the following formula: 


{ SOc =P™"( 


5 
soc = N™ sac, oO 
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Where: 


e Nbat is the number of storage units. 
e SOCo is the initial state of charge. 


As the charge and discharging operations 
are distinct, Pbess(t) can be expressed as 
follows: 


PP(t) = P(t) - BP™(t) (6) 
where: 
*a is the efficacy of charging. 
* B represents the discharge efficacy. 
* Pbc(t) is the amount of electricity used to charge 
the batteries. 
¢ Pbd(t) is the power supplied by the batteries to 
meet the demand for energy. 
As charging and discharging cannot occur 
simultaneously, Pbc(t) and Pbd(t) can be expressed 
as follows: 


{, PHOQ=e OO, 


4 (2)= (1—2x(2) 1, YP, 
Where: 


e = Ibc(t) is the efficacy of charging. 

e = Ibd(t) represents the discharge efficacy. 

e Pbc(t) is the amount of electricity used to 
charge the batteries. 

e Pbd(t) is the power supplied by the 
batteries to meet the demand for energy. 


e Vb is the batteries’ voltage 


e x(t) is a binary variable taking the value 0 
during charging phase, and the value 1 
during discharging phase. 

As charging and discharging cannot occur 
simultaneously, Pbc(t) and Pbd(t) can be expressed 
as follows: 


P'<(t)P>*(t) = 0 (8) 


This equation can be rewritten using the following 
equations, which involve the binary variable x(t). 


i x b — be 
{ (t)P P**(t)>0 6 


(1— x(t))P* - > P**(z)>0 


max 


2.4. System modelling 

Renewable energies are often used in our 
everyday lives to provide electricity for residential, 
commercial, and transportation purposes, with the 


aim of fostering the development of 
environmentally sustainable urban areas. 

The present study aims to use photovoltaics and 
wind energy in order to replicate a charging station 
for electric cars, as seen in Figure 7. The direct 
current (DC) is first derived from the solar panels 
and wind turbines. The electrical energy generated 
by the solar panel was directed towards the highest 
Power Point Tracking (MPPT) system, which 
facilitated the transmission of the highest power 
output from the solar cells. The voltage and current 
were measured in order to determine the best value. 
If the voltage fluctuates, it may be adjusted using a 
DC-DC converter. The DC Microgrid serves as a 
means of transmitting energy from its source to 
electric vehicles (EVs). The DC-AC inverter is 
used at the receiving end to facilitate the conversion 
of the current. The LC filter serves the purpose of 
attenuating undesired noise present in the circuit 
and enhancing the overallwaveform integrity of the 
electrical signal. 

The wind turbine is linked to DC-DC converters, 
whereby the electric current is monitored using 
Maximum Power Point Tracking (MPPT) 
technology in conjunction with a charge controller. 
Simultaneously, the battery bank underwent 
charging and discharging processes inside the DC 
microgrid via the use of a bidirectional DC-DC 
converter. This article mostly examines solar 
panels. 


DC Microgrid 


Bi-Directional 
DC/DC Converter 


Charge 
Controller 


Figure 1: Diagram for charging station using 
renewable energy 


2.5. Load Profile 


The electrical load consists to charge the full 
batteries of 2 cars per day is 121 kWh per day: 
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3. PROBLEM FORMULATION 


Based on the functional constraints, we present the 
formulation of the constrained sizing problem in 
this section. The following section describes the 
objective function and the various constraints. 


3.1. Objectif function 
The objective function represents the overall cost 
of installing the system, denoted as Ct. 


C.= Y)_,((N*)(h + Ry + OMY) (10) 


Where 


k refers to the system’s components, such 
as: Component | is the PV _ panel, 
Component 2 is the Wind micro turbine 
and Component 3 are the batteries 

Ik is the Initial cost of component k. 

Rk is the annual replacement cost of 
component k. 

OMk represents the yearly expenses for 
operating and maintaining component k, 
amounting to 1.5% - 2% of the initial 
investment. 


3.2. Constraints in continuous form 
The optimization problem is constrained by the 
subsequent factors: 


1. Demand constraint refers to the 
requirement that the energy generated 
from the power system components 
must be sufficient to satisfy the 
current demand. 


D(t)= PP(t)+ P¥(t)- Poe(t)}+ Pht) (11) 


Where: 
e 


D(t) is the demand in function of time. 


Ps (t) is the forecasted PV panel power 
production in function of time. 


pwf (t) is the forecasted wind power 
production in function of time. 


Npv is the number of PV modules 
Nw is the number of wind turbines 


2. State Of Charge Bounds: To safeguard 
the batteries against excessive and 
insufficient charging, adherence to the 
subsequent restriction regarding the 
highest and lowest possible values 
(SOCmax and SOCmin) is required. 


Table 1 : Net energy needed and loading time for both 


EV models 
Daily Net 
Electric Energy 
Vehicle Net energy Daily Loading Needed 
Model Needed (kWh) | Time (hour) (kWh) 
Renault 
Zoe 22 3 66 
Hyundai 
Kona 11 5 55 
Total Daily 
Total EVCS net energy 
power: 33 Average loading needed: 121 
kWh time: 4h kWh 
Nb SOC™" < SOC(t)< N"SOC™* (12) 
With: 
e Ssoc™™* is the maximal the Maximum 


leg —x(t))P* na. > P**(#)> 0 


Where: 


Recommended Operating Current (SO for 
the storage unit. 


SOC™" is the minimum amount of energy 
that ought to be retained in the storage 
unit. 

3. Ramp constraint: Rapid battery 
charging and discharging can have 
detrimental effects on the system's 
performance and shorten the batteries’ 
lifespan. As a result, the charge and 


discharging capabilities are 
constrained by the subsequent 
inequalities: 
a ib: 
r> || and r>| | a3) 
ad dt 


ris the maximal ramp rate. 

4. Complementary constraint: The 
expression for the complementarity 
between the charging and discharging 
operations is 


z(t) P? na > P(t) > 0 i 


PP ax is the maximal exchanged power by 
the storage system 

x(t) is a binary function such as x(t)=1 
during charging process, and x(t)=0 during 
discharging process. 


eee 
3404 


Journal of Theoretical and Applied Information Technology 
30" April 2024. Vol.102. No 8 
© Little Lion Scientific 


.. =m 
a! 


ISSN: 1992-8645 www. jatit.org E-ISSN: 1817-3195 


5. Space constraint: Given the following NPY as 2 Ni > 0 
inequalities, the space allocated for Net >N.> 0 
' : max — ¥2=— 
the system's components is ae 
constrained. Ni" wax 2 Nz = 0 


N?* ax = Ni> 0 NOY PT ( — BP Pe > Dy 


NY 2 N2> 0 (15) N™t(5) —Siin) > y(eere — ad P*) & > N™*(S)—S.ax); k = 1...N 
neat 2 Nz > O = 

- x Pt = PM > 0 

Where: (1 — xP ee 2 PS 

rit > P*— P, 3°* > — rit 


be be 
e WP. is the maximum quantity of rét = Py — Pia’ = — rot 


ermitted PV panels. 
p p (18) 


e W*®',,.. is the maximum quantity of wind 
turbines permitted. 4. OPTIMIZED MODELLING 


e jt__. represents the maximum quantity 

of permitted batteries. The — optimization — problem 
formulation generates a MILP problem, 
which is an NP-hard problem that cannot 
be solved by a polynomial-time algorithm. 
The resolution time of such problems is 
predominantly determined by the number 


A. Functional constraints expressed as 
discrete variables 
To implement simulation of our model, it is 


imperative to employ a discrete form that of binary variables included in the 

handles a restricted set of values. These values formulation, the geometry of the problem, 

correspond to a finite number of iterations, and the number of variables. 

referred to as sampling times t;, which are In order to mitigate potential 

defined as follows: challenges related to resolution times and 
eid tf. =a i0...N (16) accuracy arising from the expanding scale 


Ww of the MILP problem, a multitude of 
resolution algorithms have been devised. 
These include the Branch and Bound, 


© to is the initial time Genetic Algorithm, Ant Colony 


© — tris the final time Algorithm, and Particle =Swarm 

@ Nis the horizon of the study. Optimization. A summary of a few of 
The sampling interval is defined by: these techniques is presented in [18]. 

The two algorithms utilized in this 

tie ti _tF to (17) paper are the Branch and Cut Algorithm 

, 2 a) and the Genetic Algorithm. The following 

The process of discretizing the equations [10-15] software applications may be utilized to 

produce the subsequent expressions: implement these algorithms: Matlab, 


Cplex, Tomlab, Gurobi, and _ Knitro. 
Certain software applications may even 
integrate multiple algorithms during the 
resolution procedure. 

We propose a solution to this 
issue utilizing the Genetic Algorithm and 
the Branch and Cut Algorithm in the 
following section. 
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5. CASE STUDY py OW 
Prey 1000W 
To evaluate the efficacy of the Pina ay 
algorithms in addressing the sizing issue and So 4000W/h 
determine the most optimal solution, we ™ 400W/h 
conducted comprehensive numerical — 4500W/h 
simulations, which will be elaborated upon in Te, Td 3000W/h 
the following sections. The simulations were ae ! 
executed via the Tomlab user interface. ata 1,2 


According to the climatic data Table 3: Products With Their Corresponding Prices 


gathered in the city of Tangier in 2022, 
December experiences the most severe climatic 


Product Price 


conditions. As a result, we incorporated CSUN 255-265-60P | 158$ 
December's meteorological data, including | 12V 265WP PV panels 
solar radiation, ambient temperature, and wind Micro wind turbines | 285$ 
speed, into our simulations. Furthermore, the ya SA 12V/24V 


electrical power consumption values of the 
street lighting luminaires that were documented 
during the identical time were incorporated. 
Each of the datasets comprises four values 
per day: 7:00 am, 01:00 pm, 7:00 pm, and 


Vrla sealed lead acid | 316$ 
batteries 12V 230Ah 


3.3. Results 


01:00 am. The first time is at 1:00 pm on 
December 1, 2022, and the last time is at 1:00 
EF _fo . : 
yw IS the sampling 
interval, and N = 124 is the horizon. 
The values utilized in Table I were 
incorporated into the simulations; on the other 


pm on January 1, 2023. 


The subsequent are the outcomes derived 
from the simulation that incorporated the 
Branch and Cut algorithms alongside the 
Genetic Algorithm: 
The results suggest that the Genetic Algorithm 
requires a considerably longer period of time to 
converge than the Branch and Cut Algorithm. 


hand, Table II details the production and Moreover, it is evident that the acquired 
storage units implemented in the simulations. solution lacks optimality and __ further 
(iI). refinement is required. 

The methodology employed, along with the 
advantages and disadvantages of the present 
strategy, are detailed in Table 4. Eight 


Table 2: Parameters Values For The Simulation 


Parameter Value scholarly articles were subjected to our 
NPY max 15 Units examination, each employing unique research 
N™ max 6 Units methodologies across multiple levels. 

Nb* ax 5 Units 


Table 4: The Metrics And Drawbacks Of Alternative Methodologies Employed In Prior Investigations 


Author/ Methodology Advantage 
Petrusic et al. / Multicriteria | e Maximized the use of 
optimization [30] renewable sources. 
e Flexibility in charging. 
e Easy to implement. 
Noman et al. / Interval-based | e Cost Effective method. ° 
approach [31] e Reduced greenhouse gas factors. 
emissions and promote |e Required significant investment in 
sustainable transportation. wind turbines and other 
e Provided a constant power infrastructure. 
supply for specific time 
intervals. 
e Suitable for urban areas with 
high demand for EV charging. 


Disadvantage 

e Required significant investment in 
infrastructure, such as solar power 
plants and energy storage systems. 


Feasibility depends on_ several 
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Sakib et al. / Matlab/Simulink [32] 


Reduced the negative impact 
on the environment and 
geological aspects. 

Meet the demand for the next- 
generation power system. 
Better design and 
optimization of the system. 
Enhancing fuel efficiency. 


Required initial investment and 
infrastructure development. 
Depends on the availability and 
accessibility of natural gas 
resources. 


Foley et al. /PLEXOS [33] 


Minimize the cost of 
electricity dispatch and to 
some constraints. 

More accurate analysis of the 
impact of electric vehicle 
charging under peak and off- 
peak charging scenarios. 


May not be directly applicable for 
other electricity markets. 


Abas et al. / Microbial Fuel Cell 
[34] 


Developed more efficient and 
stable water-splitting systems. 
A promising technology for 
generating clean and 
renewable energy. 

Can withstand harsh 
conditions and resist 
corrosion. 


Used expensive catalysts. 

Need for large areas of land and 
the requirement for continuous 
monitoring and maintenance. 


Sayed et al. /DC microgrid [35] 


Aboelezz et al. / Generic electric 
vehicle battery charging system and 
ANN [36] 


Reduced carbon emission 
Provide sufficient electricity 
supply. 
Accurate and efficient model. 
Achieved the maximum 
generated power from the 
wind turbine across the day. 


Costly and may _ require a 
significant initial investment. 


Required a full prediction of the 
wind signature. 

Required initial investment in wind 
turbines and charging 
infrastructure. 


Karmaker et al. Fuzzy inference 
system [37] 


Reduced energy costs 
Provided lower charging costs 
Establish a sustainable 
charging infrastructure 
Effective energy management 
system. 


Effectiveness may vary with sizes 
and locations. 
Required initial investment costs. 


6. CONCLUSION 


Population growth, greenhouse gas 
emissions, and depletion of nonrenewable 
resources have all contributed to the emergence 
of novel research directions in renewable 
resources. This study seeks to comprehend the 
applications of emergent technologies for 
extracting renewable energy. In order to reduce 
greenhouse gas emissions through the efficient 
utilization of renewable resources, the 
proposed study is beneficial. The main 
emergent technologies are comprised of 
numerous sub-technologies. To purify the 
planet, it is essential to support these renewable 
technologies. Furthermore, it is worth noting 
that renewable energy sources are gradually 
surpassing the cost-effectiveness and efficiency 
of non-renewable energy sources, which have 
historically demanded substantial investments 
and operated at a relatively high level of 


efficiency. The current worldwide transition 
towards renewable and sustainable energy 
sources, as well as their implementation in 
electric vehicle charging stations, is motivated 
by cost reductions, environmental 
considerations, and the pursuit of enduring 
sustainability. 

As a result of Morocco's substantial energy 
consumption, electric vehicle (EV) recharge 
stations may be equipped with these renewable 
energies. 


A methodology is outlined in this paper 
for determining the optimal dimensions of a 
self-contained hybrid system that is designed 
for installation at an electric vehicle recharge 
station in Tangier. The sizing problem 
formulation was executed through the 
utilization of the MILP approach. 
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